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Abstract 

We present numerical simulations of the dynamics of two-dimensional Josephson 
junction arrays to study the mechanism of mutual phase locking. We show that in 
the presence of an external magnetic field two mechanisms are playing a role in 
phase locking: feedback through the external load and internal coupling between 
rows due to microwave currents induced by the field. We have found the parameter 
values (junction capacitance, cell loop inductance, impedance of the external load) 
for which the interplay of both these mechanisms leads to the in-phase solution. 
The case of unshunted arrays is discussed as well. 
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Josephson junctions are natural voltage-controlled oscillators [1]. By building up one- 
dimensional (ID) and two-dimensional (2D) arrays of the junctions one can hope to 
obtain the necessary output power which is required for many practical applications. 
This goal will be achieved when all junctions of the array are phase-locked. While a 
detailed theoretical analysis of phase locking in ID systems has been performed several 
years ago [2], [3] and has been supported by promising recent experiments [4] the full 
examination of the problems connected with the formation of stable phase locking in 2D 
arrays has not been achieved up to now. Several authors have shown that 2D arrays are 
more stable against non- uniformities in the critical currents [5]- [8]. Nevertheless, it is not 
understood so far whether 2D arrays possess an internal mechanism of coupling which 
could be a potential advantage. Even a qualitative observation regarding the mechanism 
could provide useful criteria for the problem of optimum array design. 

In this letter we address these issues by means of simulations of the dynamical proper- 
ties of 2D Josephson networks with and without external load in presence of an external 
magnetic field. 

We investigate arrays consisting of active junctions connected by superconducting 
inductances in the direction perpendicular to the bias current. Such a structure has been 
accepted as a standard design now [8]- [10]. We have N rows of M junctions. The behavior 
of every Josephson junction is simulated using the RSJ model. Within this model the 
current through the kth junction is given by 

PcVk + ¥k + siny2 fe = i k , (1) 
i k = h/Ic- 

Here (3 C = 2%CIcRn 2 /*&o is the McCumber parameter, ip k is the Josephson phase 
difference over the junction, $ stands for the magnetic flux quantum. C, Rn and Iq are 
the junction capacitance, normal resistance and critical current, respectively. I k is the 
total current flowing through the kth junction. Introducing dimensionless parameters, we 
can formulate flux quantization conditions for every cell containing junctions m and n: 

Vm - Vn = Ve J m > n) + /oE (2) 

P =i 

X> P = 1, (3) 
P =i 

where we introduced the parameters l = 27rI c L /Q and tpex^" 1 ^ = 27r$ ea; ^ m ' n ^/$o- 
Furthermore, i^ 1 '^ = I^ m ' n ^/I c is the normalized current flowing through the respective 
inductance of the pth branch. The quantity ji p is the contribution to the loop inductance 
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l from the pth branch of the cell. The frequency of the Josephson oscillations is 

Uj = Vj] (4) 

here Vj indicates the normalized DC component of the voltage across junction j. We 
restrict our simulations to the case of a uniform external magnetic flux op ext ^ m ' n ^ — <Pext- 
The resistor in the external load r s has been choosen to match the array normal resistance 
r array = Nr N /M. The first and second harmonics of the current flowing through the 
external load have been calculated as 

inu — 1/T / i(t) exp (— incut) dt, n—1,2 (5) 
Jo 

with the averaging time T > 1000/cj. In the simulations we also introduced a small spread 
of the critical currents a of about 5 % with the product V c = I c Rn beeing the same for all 
junctions. The arrays have been investigated by means of the PSCAN program [11], [12]. 

The most effective phase locking mechanism for both ID and 2D arrays is high- 
frequency coupling between an array and external load, as has been discussed by sev- 
eral authors before [1]. In dependence on inductive or capacitive character of the load 
one can expect different types of coherent solutions. The in-phase state of a ID array is 
characterized by conditions 

ifi -(pj = (6) 

for the phase difference between any two junctions. This and other types of coherent 
solutions for the ID have been described in the papers [3], [13]. For the 2D array with 
ip ext = and all junctions beeing identical we can write down the same relations for any 
phases in both "transverse" (first index) and " longitudial" (second index) directions 

<Pi,k ~ <Pu = 0. (7) 



The situation changes when we consider 2D structures with nonzero magnetic field 
within the cells. In this case the array symmetry is broken and we can no longer expect 
uniform oscillations of all junctions within a row. However, we still can define an in-phase 
solution 

¥i,k ~ = 0. (8) 

with respect to the "longitudial" direction. In the following we have to find conditions 
under which an external magnetic flux does not disturb this in-phase solution. We show 
that the 2D array standard design itself can not assure this solution and must be improved. 
For the sake of clarity we start with the four junction array (Fig. 1) and the traditional 
design of the unit cell. In this case the additional branch A is absent. Only the external 
load B must control the phase locking of the whole system. 

In Fig. 2 we show the results of simulations which have been performed for the 
value l = 0.5 of the loop inductance and two different types of external load. For the 
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capacitive load (solid line in Fig. 2a) we have obtained the anti-phase solution within 
the whole range of the external magnetic field. We find no first harmonic in the net 
oscillation. Consequently, the microwave current through the load is characterized by 
the second harmonic component essentially (Fig. 2b). The phase relations between the 
junctions in each column remain the same for all values of the external field. On the 
other hand, in the inductive load regime (filled triangles in Fig. 2a) we have got in-phase 
oscillations within the regions ip ex t < 2.0 and tp ext > 4.2 and the anti-phase solution for 
2.0 < ip ext < 4.2. It is worth noting that an increase of the magnetic flux (0 < <p ext < ir) 
leads to an increasing oscillation frequency. As a result, the transition observed can not 
be caused by the action of the external load which remains in the inductive regime. 

In order to explain the result we turn to the properties of the unshunted array first. 
Here, the microwave currents induced by external magnetic fields can contribute not only 
to the phase locking within one row but also to the phase locking between neighbouring 
rows. A stability analysis for zero Josephson junction capacitance f3 c = and Iq <C 1 has 
been performed analytically [14]. As a result one observes that even small magnetic fields 
favour the stability of the anti-phase state. 

For non-zero capacitance of the Josephson junctions and for the unit cell inductance 
Z ~ 1 the frequency dependence of both in- and anti-phase states is more complex. The 
boundaries of the stability regions can be determined numerically only. In Fig. 3 we show 
these regions for the network without external load and ip ext = 1.5. One recovers that the 
point l = 0.5 and f3 c = 0.0 (which are our parameters for Fig. 2) lies within the region 
where the unshunted array exhibits a stable anti-phase solution. 

The state choosen by the array with shunt depends on the interplay of the two high- 
frequency currents, flowing through a given junction of the array, i.e. the shunt current 
and the circulating current. For a small value of the external field the shunt current is 
strong enough to support the in-phase state. The increasing magnetic field within the 
cells produces a shift of the oscillation phases along the vertical direction and leads to a 
decreasing of the microwave current through the external load (Fig. 2b). The mechanism 
providing the in-phase solution weakens and for some value ip 1 ^ fails. 

Our simulations for the shunted array with a larger number of junctions (up to 6 in 
every row) have shown the following: The value of the critical flux (p 1 ^ for which switch- 
ing into anti-phase state takes place becomes smaller with the number M of Josephson 
junctions per row growing. From a practical point of view this means that the increasing 
of the array size makes the in-phase state less tolerant toward the external field. For a 4 
x 4 shunted array with the loop inductance l — 1.0 this limiting value is ip 1 ^ = 0.7. Af- 
ter switching the array reveals a different behavior: The phase differences of the junction 
voltage oscillations are equal to n for nearest neighbours in each column. Large microwave 
currents are flowing along the transverse inductances. 

Another conclusion we have made from the simulations of the arrays with M, N > 2 is 
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that the boundaries of stability of the in-phase and anti-phase solutions for the unshunted 
array at the /o,/3 c -plane do not depend on M and N. One possible explanation is that 
the circulating microwave currents choose the nearest ways to flow and consequently only 
the short-range coupling between neighbouring cells determines the stability conditions. 

The mechanism assuring the in-phase state depends mainly on the impedance in the 
vertical direction (in contrast to the other one created by circulating currents). This 
way one can hope to influence the relative strenght of both mechanisms by changing the 
impedance of the cell into the direction perpendicular to the bias current. We have done 
so by including the additional branch A. Fig. 4 presents the boundaries of the in-phase 
and anti-phase solutions for the unshunted (stars) and shunted (filled boxes) arrays in 
comparison. The additional branch A improves the stability of the in-phase state within 
the region l < 1.0, f3 c < 1.0 for the unshunted array (this set of parameters is supposed 
to be the most promising for obtaining maximum output power). As a result, the shunted 
array remains phase locked for all values of the external flux (Fig. 5). 

In conclusion, we have shown that the magnetic field in two-dimensional arrays pro- 
vides an internal mechanism of phase locking. This mechanism yields the stability of 
the anti-phase state for the array with inductances only beeing used for the transverse 
connections. A capacitive shunt parallel to these inductances improves the stability of 
the in-phase state. 

The authors would like to express their thanks to Deutsche Forschungsgemeinschaft 
for supporting this work. 
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FIGURE CAPTIONS 



Fig. 1. An array of 2x2 Josephson junctions with an additional branch r x , c x (denoted 
by A) and external shunt r s ,c s ,l s (B); l x = (l/4)/ , l y = (1/4)Z - 



Fig. 2. (a) The frequency w of a shunted array with inductive (filled triangles) 
and capacitive (solid line) shunt as a function of the external flux (p ext (Parameters: 
i dc = 1.5, P c = 0.0, l = 0.5, l s = 1.0, r s = 1.0). (b) The first and second harmonics of 
the shunt current i u (solid line), ( dot line) for the inductive shunt (c s = 2.0) and the 
second harmonic (open triangles) of the shunt current for the capacitive shunt (c s = 0.2) 
as a function of the external flux ip ext (Parameters: idc = 1-5, (3 C = 0.0, Zo — 0.5, l s = 
1.0, r s = 1.0). 



Fig. 3. The boundary of the in-phase and anti-phase solutions at the (3 C , / -plane for 
the unshunted array from Fig. 1 (Parameters: id c = 1-5, (p ex t = 1.5). 



Fig. 4. The boundary of the in-phase and anti-phase solutions at the (3 C , / -plane for 
unshunted array (stars) and shunted array (filled boxes) with the additional branch A 
(Parameters: id c = 1.5, ip ext = 1.5, l s = 1.0, r s = 1.0, c s = 2.0, c x = 6.0, r x = 0.1). 

Fig. 5. The frequency uu of the shunted array without (stars) and with (solid line) 
additional branch A as a function of the external flux tp ext (Parameters: idc = 1-5, f3 c = 
0.0, l = 0.5, l s = 1.0, r s = 1.0). 
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